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ABSTRACT 



Modern digital crosscorrelators permit the simultaneous measurement of all 
four Stokes parameters. However, the results must be calibrated to correct for the 
^ I polarization transfer function of the receiving system. The transfer function for any 

ly-^ I device can be expressed by its Mueller matrix. We express the matrix elements in 

f^ ' terms of fundamental system parameters that describe the voltage transfer functions 

f"^ ■ (known as the Jones matrix) of the various system devices in physical terms and 

thus provide a means for comparing with engineering calculations and investigating 
the effects of design changes. We describe how to determine these parameters with 
O,. astronomical observations. We illustrate the method by applying it to some of the 

JL ! receivers at the Arecibo Observatory. 

J^ ' Subject headings: polarization — instrumentation: polarimeters — techniques: 



polarimetric 



1. INTRODUCTION 

The polarization response of a single-dish radio telescope is usually described in terms of the 
native polarization of the feed. For example, a dual-linear polarized feed provides two outputs 
("channels") which are orthogonal linear polarizations; the sum is Stokes / and the difference 
Stokes Q. Correlating the two outputs with zero phase difference provides the other linearly 
polarized Stokes parameter U, and correlating with 90° phase difference provides the circularly 
polarized Stokes parameter V. Similarly, with a dual-circular feed the sum and difference provide 
/ and V, and the correlated outputs provide Q and U. However, these statements are only 
approximate. Feeds are almost never perfect, so their polarizations are only approximately 
linear or circular; moreover, some feeds, such as turnstile junctions, can have two orthogonal 
polarizations with arbitrary ellipticity that is frequency-dependent. And generally speaking, no 
feed has two outputs that are perfectly orthogonal. 



The feed's response is modified by the electronics system, which introduces its own relative 
gain and phase differences between the two channels. These must be calibrated relatively 
frequently because they can change with time. For example, if the feed is perfect and the relative 
gain differs, then the difference between the two channels is nonzero for an unpolarized source, 
making the source appear to be polarized. For a perfect linearly polarized feed, if the relative phase 
of the two channels differs, then a linearly polarized source appears to be circularly polarized. 
Thus two quantities, the relative gain and phase, must be calibrated. This is most effectively done 
by injection of a correlated noise source into the two feed outputs or, alternatively, by radiating a 
noise source into the feed with a polarization that provides equal amplitudes in the feed outputs. 

The modification of the four Stokes parameters by these system components is most generally 
described by a 4 x 4 matrix. This matrix is known as the Mueller matrix (Tinbergen 1996). Along 
with every Mueller matrix goes a Jones matrix, which describes the transfer function for the 
voltages. We express the Jones matrix elements in terms of complex voltage coupling coefficients. 
We consider three such matrices: one describes the polarization state of the feed (circular, 
elliptical, or linear); one the nonorthogonality of the polarized outputs; and one the relative gain 
and phase of the electronics system. We calculate the associated Mueller matrices and describe 
how to solve for the coupling coefficients from astronomical measurements. We illustrate the 
method by applying it to some of the systems at Arecibo Observatory^. 

We begin by reviewing the basic theoretical concepts and developing the structure of our 
treatment. §g introduces the Stokes parameters, Mueller matrix, and Jones vector and matrix. §^ 
defines the Mueller matrices for the different effects we describe, including the mechanical rotation 
of the feed with respect to the sky, the imperfections in the feed, the feed, and the amplifier chain. 
§Q discusses the combined effects of the receiver components and calculates that matrix product. 
§^ describes the technique for evaluating the matrix elements from observations. §^ describes how 
to apply the derived Mueller matrix to correct for the instrumental effects. Finally, §|^ provides 
illustrative results for the two 21-cm line receivers at Arecibo Observatory. This paper is a more 
general and succinct version of Arecibo Technical and Operations Memo 2000-04, which covers all 
of Arecibo's receivers. 



2. STOKES, MUELLER, AND JONES 

The basic reference for our discussion of the fundamentals is the excellent book on astronomical 
polarization by Tinbergen (1996). A more mathematical and fundamental reference is Hamaker, 
Bregman, and Sault (1996), which the theoretically-inclined reader will find of interest. In the 
following we make several unproven statements and assertions about Mueller and Jones matrices; 



^The Arecibo Observatory is part of the National Astronomy and Ionosphere Center, which is operated by Cornell 
University under a cooperative agreement with the National Science Foundation 



the explanations and justifications can be found in the abovementioned references. 



2.1. The Stokes and Jones vectors 

The fundamental quantities are the four Stokes parameters (/, Q, U, V), which we write as the 
4-element Stokes vector 



Q 

u 

V 



The Stokes parameters are time averages of electric field products; we use the terms "voltage" and 
"electric field" interchangeably because the radio telescope's feed converts one to the other. The 
Jones vector represents the fields as orthogonal linear polarizations {Ex,Ey), and considers them 
as complex to account for their relative phase: 



Ex 

Ey 



(2) 



Instructive special cases include pure linear and pure circular polarization. Orthogonal 
linear polarizations are, obviously, (Jx, Jy) = ([1,0], [0, 1]) (In the text, we write these vectors 
as rows instead of columns for typographical purposes). Orthogonal circular polarizations are 
(Jl, Jr) = (~^) /^)' '^ti^re i = \/— 1. Orthogonal polarizations satisfy, for example, JxJy = 0, 
where the bar over a symbol indicates the complex conjugate and all products are time averages. 

It is straightforward to relate the Stokes parameters to the components of the Jones vector: 



/ ^ Ex Ex + EyEy 

Q = ExEx — EyEy 

U = ExEy + ExEy 

iV = ExEy — ExEy . 



(3a) 
(3b) 
(3c) 
(3d) 



2.2. The Mueller and Jones matrices 

When the fields pass through some device, such as a feed or an amplifier, they suffer amplitude 
and phase changes. These modify the Stokes parameters. The MueUer matrix is the transfer 
function between the input and output of the device: 



Sout — M • Sin • 



(4) 



The Mueller matrix is, in general, a 4 x 4 matrix in which all elements may be nonzero (but they 
are not all independent). In the usual way, we write 



M 



mil miQ 

mQi mQQ 

mui rnuQ 

mvi rnvQ 



miu miv 

rriQu rngv 

muu iTiuv 

mvu innyy 



(5) 



The matrix elements are just the partial derivatives, for example 



dVr 



mvQ 



out 



dQi 



-'-ini^ 'I 



(6) 



zri)'-' m 



Every Mueller matrix has its Jones matrix counterpart; the Jones matrix is the transfer 
function for the voltages. We defer further discussion of Jones matrices until our treatment of 
three specific cases of interest for radio astronomical systems. 



THE RADIOASTRONOMICAL RECEIVER COMPONENTS AND THEIR 

MATRICES 



In this section, we consider the Mueller matrices of devices that are encountered by the 
incoming radiation on its way from the sky to the correlator output. We consider three devices. 
The first device encountered by the incoming radiation is the telescope (Mgky), which mechanically 
rotates the feed with respect to the sky. The second device is the feed, which we split into two 
parts. The first part (Mf) has the ability to change the incoming linear to any degree of elliptical 
polarization; by design, feeds are intended to produce either pure linear or pure circular, but in 
practice the polarization is mixed, i.e. elliptical. The second part (Mift) describes imperfections 
in such a feed, specifically the production of nonorthogonal polarizations. The third device is the 
amplifier chain (Ma). 

We describe the Jones matrix of the incoming radiation in linear polarization, as in equation]^. 
However, the feed matrix Mf can radically change the polarization state; for example, as we shall 



see a dual-circular feed changes the order of the Stokes parameters in equation |l[ Thus the signal 
voltages, after going through Mp, are not intuitively described as {Ex,Ey), because the {X,Y) 
connote linear polarization. Therefore, for the voltages after the output of Mp, we will use the 
symbols {Ea,Eb), or simply (A,B), to emphasize the fact that the state of polarization can be 
arbitrary. 

We assume that the remaining device matrices Mipr and Ma produce closely-matched 
replicas of the input Stokes parameters because, by design (hopefully!), the imperfections in the 
feed are small. We will retain only first-order products of these imperfections. 



3.1. Mueller matrix relating the radio source to the receiver input 

Astronomical continuum sources tend to have linear polarization but very little circular 
polarization; we assume the latter to be zero. Moreover, we express the source polarization as a 
fraction of Stokes I. Thus, 



1 

^src 
^src 





(7) 



A linearly polarized astronomical source has Stokes {Qsrc, Ugrc) defined with respect to 
the north celestial pole (NCP). The source polarization is conventionally specified in terms of 
fractional polarization and position angle with respect to the NCP, measured from north to east. 
We have 



^ src — ^src COS Zi-^^lg^c 
^ src ~~ ^src Sm Zr J^gyc 

p - (o^ ^u"^ ^1/2 

J^ src — K^src ^ '-^ src) 



PAsrc = O.Stan ^ 



c/.. 



Qs 



(8a) 
(8b) 
(8c) 

(8d) 



As we track a source with an alt-az telescope, the parallactic angle PAaz of the feed rotates 
on the sky. PAaz is defined to be zero at azimuth and increase towards the east; for a source 
near zenith, PAaz ~ az, where az is the azimuth angle of the source. The Stokes parameters seen 
by the telescope are {Qsky, Usky), and are related to the source parameters by 



MsKY 



1 

COS 2PAaz 

-sm2PAaz 






sin 2PAaz 

cos 2PAaz 

1 



(9) 



The central 2x2 subniatrix is, of course, nothing but a rotation matrix. 

For an equatorially-mounted telescope, PAaz doesn't change unless the feed is mechanically 
rotated with respect to the telescope. This was the case with the NRAO 140-foot telescopeg, which 
was one of the several systems we used during the genesis of this work. This relative rotation 
between feed and telescope can produce unintended changes in the feed response to astronomical 
sources, which makes such telescopes inherently less accurate for polarization measurement. The 
140-foot telescope is the last of the great equatorial radio telescopes and with its recent closure 
a detailed discussion of these matters has become irrelevant; the remainder of the present paper 
considers alt-az telescopes exclusively. 



3.2. Mueller matrix for a perfect feed providing arbitrary elliptical polarization 

The feed modifies the incoming voltages with its Jones matrix. Suppose that the feed mixes 
incoming linear polarizations with arbitrary phase and amplitude, keeping the total power constant 
and retaining orthogonality; this makes it a perfect feed that responds to elliptical polarization. 
Following Stinebring (1982) and Conway and Kronberg (1969), we write the transfer equation for 
the feed as 



EA,out 
EB,out 



COS a 



P«X 



sma 



-^X 



sma 



cos a 



Ex,i 
Eya 



(10) 



The feed can completely alter the polarization state, so the output Jones voltages are more 
intuitively described by subscripts {A,B) instead of (X, y), which connote linear polarization. 
Here a is the amount of coupling into the orthogonal polarization and x is the phase angle of that 
coupling. For example, for a native linear feed a = and x = Oj fo^ ^ native circular feed a = 45° 
and X = 90° . Using this with equation |^, we find 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 



Mf 



1 








COS 2a sin 2a cos x sin 2a sin x 

sin 2a cos x cos^ a — sin^ a cos 2% — sin^ a sin 2x 

■ sin 2a sin x — sin^ a sin 2^ cos^ a + sin^ a cos 2% 



(11) 



Notice that in the right-bottom 3x3 submatrix, the off-diagonal transposed elements are of 
opposite sign for two of the three pairs and the same sign for one. This is not an algebraic error! 

Some instructive special cases include: 

(1) A dual linear feed: a = 0, x = 0, and Mp is diagonal. 

(2) A dual linear feed rotated 45° with respect to {X, Y): a = 45°, X = 0, and 



Mf 



10 
10 
0-100 
1 



(12) 



As expected, this interchanges Stokes Q and U, together with a sign change as befits rotation. 
(3) A dual linear feed rotated 90° with respect to {X,Y): a = 90°, X = 0, and 



Mf 



10 
0-100 
0-10 
1 



(13) 



As expected, this reverses the signs of Stokes Q and U. 
(4) A dual circular feed: a = 45°, x = 90°, and 



Mf 



10 
1 
10 
0-100 



(14) 



The combination (a = 45°, x = 90°) permutes the order of the Stokes parameters in the output 
vector, making it (/, V, U, —Q). {EaEa — EbEb) provides Stokes V , instead of the Q written 



in equation ^; in other words, it makes the feed native dual circular. With respect to linear 
polarization, a = 45° has the same effect as in case (2), namely to interchange (Q, C/) and change 
the sign of Q, because it is equivalent to a feed rotation of 45°. 

(5) A dual-elliptical feed: arbitrary a, x = 90°. If x = 90° then orthogonal linear inputs 
produce orthogonal elliptical outputs with the ellipticity voltage ratio equal to tan a (see Tinbergen 
1996, Figure 2.1). Thus a = 0° passes the linear polarizations without modification, a = 90° 
reverses the sign of the position angle, a = 45° produces orthogonal circular polarizations, and 
other values of a produce orthogonal elliptical polarizations. 



3.3. An important restriction: we set x = 90° 

Generally speaking, we prefer either pure linear or pure circular feeds, i.e. we prefer either the 
combination (a, x) = (0°, 0°) or (a, x) = (45°, 90°) in Mp. If a feed is designed to produce circular 
polarization with (a,x) = (45°, 90°) and instead produces elliptical polarization whose major axis 
is aligned with the X direction, then x = 90° but a ^ 45°. If the ellipse is not aligned with X, 
then X / 90°, but this is equivalent to having x = 90° and physically rotating the feed. Thus, 
without loss of generality, we can take x = 90° (also see §3.4). This leads to great simplification in 
Mp, whose restricted form becomes 



Mft 



10 

cos 2a sin 2a 

10 

— sin 2a cos 2a 



(15) 



3.4. An imperfect feed 

Again we follow Stinebring (1982) and Conway and Kronberg (1969), and represent the 
imperfections of a feed by the Jones matrix 



Ex,out 

Ey.out 



eie" 



£26" 



Ex,in 
EY,in 



(16) 



Here the e's represent undesirable cross coupling between the two polarizations; for example, this 
might be caused by the two linear probes not being exactly 90° apart. The (/)'s are the phase 
angles of these coupled voltages. This equation assumes that the feed is "good", meaning that we 
need retain only first-order terms in e (which makes the diagonal elements unity); however, for the 
moment we allow the phases to be arbitrary. 



After a little algebra, we find the matrix for the Imperfect Feed to be 



MiF 



1 Secos Sesin 

1 Aecos Aesin 

Secos —Aecos 1 

Sesin —Aesin 1 



(17) 



where Se cos = ei cos (pi + €2 cos (p2', Se sin = ei sin (pi + e2 sin (/)2 ; Ae cos = ei cos (pi — €2 cos (p2 ; 
Aesin = ei sin^i — e2sin(/)2. The imperfections in a good feed are completely specified by four 
independent parameters. 

The central 4-element submatrix is a rotation matrix that represents an error in position 
angle of linear polarization (its diagonal elements are unity because of our first-order expansion 
in e); in-phase mutual voltage coupling between the two probes causes an apparent rotation. Its 
off-diagonal elements Aecos are impossible to measure without calibration sources whose position 
angles are accurately known. Moreover, a small rotation can also occur because of mechanical 
inaccuracy in mounting the feed. In practice, these problems make it impossible to separate this 
factor, so one might as well assume it is equal to zero and make appropriate adjustments to the 
position angle ex post facto. 

Assuming Ae cos = is consistent with a more stringent assumption, namely that e = ei = e2 
and (p = (pi = (p2 so that, also, we have Aesin = 0. Physically, this is equivalent to assuming 
that the coupling in a linear feed arises from the two probes being not quite orthogonal and that 
the coupling in each has the same relative phase. In other words, it makes the correlated output 
non-orthogonal. 

This assumption might seem to be too restrictive because, by requiring (pi = (p2, it eliminates 
the possibility of imperfections inducing a change in the ellipticity of the polarization. However, it 
leads to no loss in generality of our treatment, because the out-of-phase coupling case is included in 
the feed matrix Mp. Our restricted case of an imperfect feed is described by just two parameters: 



MiFr 



1 


2e 


COS(p 


2esin 


1 










2ecos(^ 




1 





2e sin (p 







1 



(18) 



3.5. The amplifiers 

The two polarization channels go through different amplifier chains. Suppose these have 
voltage gain {gA,gB), power gain {Ga,Gb) = {9a^9b)^ ^^'^ phase delays {iPa,'>Pb)- The Jones 
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matrix is 



EA,out 
EbmuI 



9Ae 




iV'A 



gse 



itpl 



EA,i 
Eba 



(19) 



In practice, the amplifier gains and phases are cahbrated with a correlated noise source (the 
"cal"). Thus, our amplifier gains {Ga,Gb) have nothing to do with the actual amplifier gains. 
Rather, they represent the gains as calibrated by specified cal intensities, one for each channel. If 
the sum of the specified cal intensities is perfectly correct, then the absolute intensity calibration 
of the instrument is correct for an unpolarized source (Stokes / is correctly measured in absolute 
units). In our treatment, Ga + Gb = 2 by necessity because we deal with fractional polarizations. 

If the ratio of the cal intensities is correct, then the difference between the two polarization 
channels is zero for an unpolarized source. This happy circumstance does generally not 
obtain. However, the relative cal intensities are known fairly well, which allows us to assume 
AG = Ga — Gb 'C 1 and to carry gAdB to first order only, meaning we take gAdB = 1- With this 
first-order approximation, we have 



Ma 



1 

AG 

2 


2 

1 
















cos^ 


— sin^ 








sin'0 


cos^ 



(20) 



The incorrect relative cal amplitudes produce coupling of Stokes / into Q through the nonzero 
mQi] the difference between the relative cal and sky phases produces a transfer of power between 
the two correlated outputs, as we now discuss. 

The difference between the amplifier phases is also referred to the cal. Thus a phase difference 
V' = V'A — "0-8 represents the phase difference that exists between a linearly polarized astronomical 
source and the cal and has nothing to do with the amplifier chains. The behavior of this phase 
difference depends on the native polarization of the feed. 

For a perfect native dual linear feed, the phase of a linearly polarized astronomical source, 
modulo 180°, is independent of PAaz because the linearly polarized dipoles have no relative 
phase difference (x = 0)j the phase changes by 180° when the measured U changes sign as PAaz 
changes. Thus, ip modulo 180° is independent of PAaz- The desirable case '0 = means that 
{EaEb + EbEa) contains pure Stokes U and {EaEb — EbEa) pure Stokes V. 

For a perfect native dual circular feed, the phase of a linearly polarized astronomical source 
(and therefore tp) rotates as ±2PAaz (see section 5.2 and equation pOb ). At PAaz = 0, the condition 
i/j = produces the correctly defined Stokes parameters in equation |l^, {EaEb + EbEa) = Ugky 
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and {EaEb - 
PA 1-^ 



EbEa) = —Qsky If V' 7^ 0) then the correctly defined Stokes parameters occur at 



3.6. The correlator outputs 

Our measured quantities are four time averaged voltage products from the digital correlator: 
from autocorrelation, {EaEa, EbEb); from crosscorrelation, (EaEb, EbEa)- In these products 
we consider the second quantity to be delayed relative to the first. Each correlation function has 
N channels of delay. We Fourier transform (FT) these quantities to obtain spectra. 

The autocorrelation functions are symmetric and thus their FT's are real, with no imaginary 
components. We combine the two measured crosscorrelation functions into a single one with 2N 
channels; it has both negative and positive delays and is generally not symmetric, so its FT is 
complex. We combine these Fourier transforms as in equation |3f 



" APB 




AMD 




AB 




BA 





FT{EaEa) + FT{EbEb) 

FT{EaEa) - FT{EbEb) 

2 Re[FT{EAEB)] 

2 lm[FT{EAEB)] 



(21) 



If we have an ideal native dual linear feed and a perfect receiver, then {APB , AMB ^ AB ^ BA) 
(/, g, [/, V)sky] for native circular, {APB, AMB, AB, BA) = (I, V, U, -Q)sky 



4. THE SINGLE MATRIX FOR THE RADIOASTRONOMICAL RECEIVER 

4.1. The general case with x = 90° 

The observing system consists of several distinct elements, each with its own Mueller matrix. 
The matrix for the whole system is the product of all of them. Matrices are not commutative, so 
we must be careful with the order of multiplication. 

We express the Jones vector of the incoming radiation in linear polarization. The radiation 
first encounters the feed, producing Stokes parameters as specified by Mpr in equation |l^. Next it 
suffers the restricted set of imperfections associated with Mipr (equation |^) . Finally it proceeds 
through the amplifier chains, undergoing Ma (equation pO|) . The product of these matrices, in this 
order (Mtot = Ma • Mift • Mpr), produces the vector that the correlator sees, which we denote 
by COR. In calculating COR, we ignore second order terms in the imperfection amplitudes 
(e, AG) but, of course, retain all orders in their their phases (</), ip) and also in the feed parameter 
a. This gives 
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M 



TOT 



1 

AG 

2 

2ecos{(f) + -0) 
2esm((^ + ip) 



-2esin(/)sm2a H — ^cos2q) 2ecos( 



cos 2a 

sin 2a sin ip 

- sin 2a cos ip 





cosip 
sin-;/' 



(2e sin cf) cos 2a H — ^ sin 2a) 

sin 2a 

— cos 2a sin ip 

cos 2a cos ip 



.(22) 



The terms in the top row make / 7^ 1 for a polarized source. If one derives fractional polarization, 
for example y, then it will be in error by amounts comparable to [(e, AG) x [Q, U, V)]. For the 
weakly polarized sources we use as calibrators, these products are second order and therefore are 
of no concern. 

However, for a strongly polarized source such as a pulsar, these terms are first order. This 
can be particularly serious for timing, because polarization variations across the pulse will produce 
errors in the pulse shape. These effects can be eliminated by correcting for the Mueller matrix. 



4.2. Two important special cases 



Commonly, feeds are intended to be either pure linear or circular. For these two important 
cases we have the following, expanded to first order: 



(1) A dual-linear feed with a slight elliptical component, meaning that [a 
with 5a <^ 1. We ignore second order products involving 5a: 



(0°,90°) + (5a] 



M 



TOT.lin 



1 

AG 

2 



±^ 2ecos 



±1 



2€ cos{(p + ip) ziz25asm'ip 
2€sm{(p + 7p) ^25acosip 



2esm(p 
±25a 

cos ip =F sin ip 
sin ip lb cos -0 



(23) 



For terms with two signs, the top sign is for the 0° case and the bottom for the 90° case. 
(2) As above, but for a dual-circular feed with [a = (45°, 135°) -|- 5a]: 



M. 



TOT,circ 



1 

AG 
2 

2ecos(0 -|- ip) 
2esm{(j) + Ip) 



=F2e sin (p 2e cos 4 
T2<5a 

lb sin Ip cos Ip 
=F cos Ip sin Ip 



I AG 
^ 2 

±1 
ziz25asmip 
^25a cos ip 



(24) 



Recall that in the circular case, the order of the Stokes parameters in the output vector is 
permuted: [/, V, U, —Q]- As with the imperfect linear, the imperfections in Mift produce coupling 
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between Stokes / and {V,U,—Q), represented by the nonzero elements in the left column. The 
order of the terms in the column is independent of the feed polarization, but the Stokes parameters 
are not, so the imperfections produce different effects for the two types of feed. 



5. EVALUATING THE PARAMETERS IN THE MATRIX 



We evaluate the parameters in equation 22 using observations of a polarized source tracked 
over a wide range of position angle PAaz- The source is described by Sgrc and the Mueller matrix 



for the radiation entering the feed by Mgky, both described in section 3.1. The full Mueller matrix 
is Mtot ■ Mgky The product of this matrix with Sgrc results in a set of of four equations, one for 
each element of the observed COR vector. Recalling that we define the source Stokes parameters 
as fractional (so that Igrc = 1) and that we assume Vsrc = 0, they are expressed by the four 
equations embodied in 



APB 




1 


AKIB 
AB 


= Mtot • Msky • 


^src 
Usrc 


BA 








(25) 



In practice, we cannot reliably measure the PAaz dependence of APB. APB is approximately 
equal to the total intensity, and it is rendered inaccurate by small gain errors. Thus, in practice 
we use fractional correlator outputs. We define 



COR' 



COR 

APB 



and consider only the last three equations in equation 25: 



(26) 



AMB' 
AB' 
BA' 



Mtot • Msky 



WsTC 

^ src 





(27) 



Note that the division by APB produces errors in the other elements of COR', but these errors 
are second order because they are products of AG and/or e with quantities such as Qsrc that are 
already first order. Our whole treatment neglects second order products, so we can neglect these 
errors. 



Multiplying out equation |27|, we obtain three equations of the form 
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AMB = Aamb' + Bmib' cos 2PAaz + Cmib' sin 2PA^ 



(28) 



where the coefRcients {A,B,C) are comphcated function of the 5 parameters (a, e, 0, AG, V') and, 
also, the two source Stokes parameters {Qsrc, Ugrc) which are not known ah initio, so we have 7 
unknown parameters. We have 9 measured quantities, three (A, i?, C) for each correlator output; 
these are derived from least squares fits of (AMB' , AB' , BA') to PAaz- We use nonlinear least 
squares fitting to solve for the 7 unknown parameters. In practice, we use numerical techniques 
to obtain the relevant derivatives. An alternative fitting technique that is useful when combining 
the results from Nsrc different sources into a grand average is to express the coefficients Ajij^iBi, 
etc., in terms of the (5 + 2Ng^c) unknown coefficients, and to lump all observations of all sources 
together in one grand nonlinear least squares fit. 

Nonlinear least square fitting is often plagued by multiple minima, and the present case is 
no exception when the polarization is nearly pure linear or circular. To discuss these cases we 
temporarily assume (e, AG) = 0, which makes the mathematics more transparent. 



5.1. The nearly linear case 



Equat 


ion |25| becomes 






Qout 




' AMB' 




±1 





Uout 


K. 


AB' 


= 


±.26asm.'il) 


cosV' 


Vout 




BA 




~^26a cos V' 


sinV' 



cos 2PAaz sin 2PAaz 
- sin 2PAaz cos 2PAaz 



^src 
'-^ src 



(29a) 



which, when multiplied out, becomes 



AMB' 
AB' 
BA' 



± COs2PAsrc 
COS i/' sin 2PAsrc ± 25cK sin -0 cos 2PAs 
sin Tp sin 2PAsrc T 25a cos ip cos 2PAs 



± sin2PAs7-c 

- cos '0 cos 2PAsrc i 2Sa sin ^ sin 2PAs 

- sin V' cos 2PAsrc T 2'5a cos i/) sin 2PAs 



Pa, 
Ps, 



. COS 2PAa 
, sin2PA„ 



.(29b) 



For 5a = we have 



mB' 




± cos 2PA,rc 


± sin 2PAsrc 


AB' 


= 


COS Tp sin 2PAsrc 


— cos V" COS 2PAsrc 


BA 




sin V' sin 2PAsrc 


— sin ip cos 2PAsrc 



PsTC COS ZPAnz 
-isrc Sm ZPAaz 



(29c) 



This shows that, for 5a — > 0, one cannot distinguish between the two cases (a,'0) = (0°,V'o) and 
{a,ip) = (90°, -00 + 180°), because the two bottom rows change sign for ipQ + 180°. For these two 
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cases, the signs of {Q, U) change, which is equivalent to rotating the feed by 90°. The physical 
interpretation is straightforward: a = 90° converts Ex to Ey in equation |l^, thus changing the 
sign of Q; changing ip by 180° changes the sign of U; the combination is equivalent to rotating the 
feed by 90°. 

In practice one must deal with this problem. For a conventional linear feed, loosely described 
as two E-field probes in a circular waveguide, the combination {a,ip) = (90°, 180°) is physically 
unreasonable. However, for some feeds, in particular a turnstile junction operating in linear 
polarization, either possibility can occur and one must make the appropriate choice based on 
frequencies well away from where the feed is pure linear; for a turnstile we expect ^ ~ const. 



5.2. The nearly circular case 





Equation 


|25|is 










Vout 




" MIB 




^26a 







Uout 


RS 


AB 


= 


lb sin ^ 


cosV' 




^out 




BA 




^cosip 


sin^ 


and 


becomes 













cos 2PAr, 



sin 2PAn 



sm2PAr,., cos2PAa.,, 



^src 



(30a) 



AMB 
AB 
BA 



^25a cos 2PAsrc 
sm{2PAsrc ± V') 

^ COs{2PAsrc ± V') 



^25asm.2PAsrc 

- COs{2PAsrc ± V') 

^ sin(2P^,^c ± ^) 



^src COS Ar J\(i2, 
lorn Sm ZrA.nz 



(30b) 



This shows that the angles 2PAsrc and ip are inextricably connected. We can determine only their 
sum (for a = 45°) or their difference. It also shows that the two solutions a = (45°,— 45°) are 
degenerate as 6a — > 0, because the bottom row changes sign for these two cases (thus flipping the 
derived sign of Qsrc and thereby rotating the derived PA^rc by 90°). The physical interpretation 
of this degeneracy is straightforward: for a pure circular feed, the phase of a linearly polarized 
source rotates with 2PAaz and its absolute value depends both on the system phase ip and the 
source position angle PAgrc- 

In practice, one can only deal with this problem by having additional information, namely 
knowing either ip or PAgrc- For the case of turnstile junctions, which are narrow band devices, one 
can determine ip at frequencies well away from that where the feed is pure circular and interpolate. 
Then, during the fit, this value of ip should be fixed. For a wideband circular feed, there is no 
substitute for an independent calibration of the linearly polarized position angle, either with a test 
radiator or with a source of known polarization. 
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5.3. Commentary 

(1) e is the quadrature sum of the P^az-independent portions of the two correlated outputs 
{AB, BA). This power is distributed between those outputs according to ((/> + V')- In the near-hnear 
case, ip can change by 180° by changing the choice for a, and this also produces a 180° change in 

<t>. 

(2) Consider a high-quality standard linearly polarized feed that has a correlated cal 
connected by equal-length cables. Such a feed has a ~ 0° and the equal-length cables mean that 
ip ~ 0°. However, the solution yields ip ~ 180° if the sign of AMB is incorrect, which can easily 
happen if one interchanges cables carrying the two polarizations; this is equivalent to reversing the 
handedness of PAaz and PAgj-c- 

(3) If one has a system without a correlated cal, then ip is meaningless. ^ has contributions 
at r.f. (from the difference in cable and electrical lengths to the first mixer) and i.f. (from length 
differences after the mixer). Normally the latter is likely to dominate because the cable runs from 
the feed to the control room are long. For example, at Arecibo we found -^ ~ 0.1 rad MHz~^, 
roughly constant among different systems. This corresponds to an electrical length difference of 
~ 5 m, most of which probably occurs along the pair of ~ 500 m optical fibers that carry the two 
channels from the feed to the control room. 

(4) We have adopted the following procure for phase calibration. If there is a correlated cal 
then we measure ipcai and fit it to a constant plus a slope -4^; we subtract this fit from the source 
phase and produce corrected versions of {AB, BA). Thus, the only component left in the correlated 
products is the difference between source and cal phase, which is the same as ip in the equation ^. 

If there is not a correlated cal, then we measure ^iligrc and fit it to a constant plus a slope -$-\ 
we subtract the slope but not the constant from the source phase and produce corrected versions 
of {AB^BA). While most of this slope is in the system, this procedure also subtracts away any 
intrinsic slope caused by Faraday rotation. 



6. APPLYING THE CORRECTION 

One of the major reasons to determine the Mueller matrix elements is to apply it to 
observations and obtain true Stokes parameters. There are two steps to this process. 



6.1. Applying Mtot and M^ky 

The completely general form of equation ^ uses observed voltage products instead of 
fractional ones and does not force V = Q. Thus, to derive the source Stokes parameters from the 
data, we use 
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^src 
Wsrc 
^ src 



(Mtot • Msky)" 



APB 

MIB 

AB 

BA 



(31) 



6.2. Deriving true astronomical position angles 

The position angle of the source polarization PAgrc is defined relative to the local 
idiosyncrasies. For a dual linear feed these include the angle at which feed probes happen to be 
mounted and, also, which feed probe happens to be defined as A. For a dual circular feed this 
includes the phase angle at which the correlated cal happens to be injected and the angle at which 
the feed happens to be oriented. 

Astronomers wish to express position angles in the conventional way, viz. with PAsrc 
measured relative to the North Celestial Pole. There is also the possibility of its handedness, 
but this is taken care of automatically in the fitting process if the PAaz is correctly defined. To 
satisfy the astronomers' desire, we must apply a rotation matrix Mastron; which looks like Msky 
in equation |9| with PAaz replaced by Oastron- For a linearly polarized feed the angle Oastron is 
the angle of the feed probes with respect to azimuth arm. There is a sign ambiguity that is best 
determined by empirical comparison with known astronomical position angles. At low frequencies, 
one must include the effects of terrestrial ionospheric Faraday rotation, which is time variable. 



7. SAMPLE RESULTS 

Here we present sample results for the two Arecibo L-band receivers (L-band is the frequency 
range surrounding the 21-cm line). For each observing session the digital correlator was split into 
four 25-MHz chunks centered at different frequencies. The L-band Wide receiver (LEW) is a very 
wide-bandwidth dual linearly polarized feed. The L-band Narrow receiver (LBN) is a turnstile 
junction whose polarization state changes from dual linear to dual circular over a frequency range 
~ 100 MHz. First, however, we reiterate the definitions of the parameters. 



7.1. Reiteration of parameter definitions 

AG is the error in relative intensity calibration of the two polarization channels. It results 
from an error in the relative cal values (TcaiAjTcaw)- Our expansion currently takes terms in 
AG to first order only, so if the relative cal intensities are significantly incorrect then the other 
parameters will be affected. 
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The relative cal values should be modified to make AG = 0, keeping their sum the same. To 
accomphsh this, make TcaiA,m.odified = TcaiA (l - ^) and TcaW, modified = TcaW (l + ^) • 

■j/) is the phase difference between the cal and the incoming radiation from the sky; see the 
discussion following equation ^. It redistributes power between ([/, V) for a dual linear feed and 
between (Q, U) for a dual circular feed (equations p3 and [2? 



Q is a measure of the voltage ratio of the polarization ellipse produced when the feed observes 
pure linear polarization. Generally, the electric vector traces an ellipse with time; tana is the 
ratio of major and minor axes of the voltage ellipse. Thus, tan^ a is the ratio of the powers. If a 
source having fractional linear polarization Psj-c = \/Qirc + ^'src i^ observed with a native circular 
feed that has a = ^ + Ja, with (5a <C 1, then the measured Stokes V will change with 2PAaz and 
have peak-to-peak amplitude Ada. 

X is the relative phase of the two voltages specified by a. Our analysis assumes x = 90°; this 



incurs no loss of generality, as explained in sections p^ and 3.4. 



e is a measure of imperfection of the feed in producing nonorthogonal polarizations (false 
correlations) in the two correlated outputs. Our expansion takes e to first order only. The only 
astronomical effect of nonzero e is to contaminate the polarized Stokes parameters (Q, U, V) by 
coupling Stokes / into them at level ~ 2e; the exact coupling depends on the other parameters. 
For weakly polarized sources, this produces false polarization; for strongly polarized sources such 
as pulsars, it also produces incorrect Stokes /. 

4> is the phase angle at which the voltage coupling e occurs. It works with e to couple / with 

^astron is the angle by which the derived position angles must be rotated to conform with 
the conventional astronomical definition. 



7.2. Results for LBW 

LBW is a very wide-band feed with native linear polarization. It has some problems with 
resonances. Apart from these, the parameters are nearly ideal and frequency-independent from 
1175 to 1680 MHz: a ~ 0.25° and e ~ 0.0015. Care was taken by the receiver engineer to equalize 
the length of the two cables for the correlated cal; as a result, ip ~ 4.6°. 

Figure || exhibits observational results for two sources obtained during two Mueller-matrix 
calibration observing sessions. The digital correlator was split into four 25-MIIz chunks, with 
duplicate coverage at 1415 MHz on the two days; this provides seven different center frequencies. 
The change of position angle from Faraday rotation is obvious over the broad band covered by the 
data, and even within a single spectrum. Regarding fractional polarization, the decrease toward 
low frequency for B0017+154 is not surprising; this is caused by Faraday depolarization. However, 
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that of B1634+269 reaches a niinimum near 1400 MHz, and this behavior is somewhat unusual. 
We beheve that our measurements are correct and that this observed behavior is real. 



7.3. Results for LBN 

LBN is a turnstile system without a correlated cal and is commonly used over a large frequency 
range. Turnstile junctions are narrow-band devices for two reasons: one, the polarization response 
is defined by physical path length differences in waveguide, and these length differences correspond 
to incorrect phase differences as one moves away from the design frequency; two, unwanted 
reflections within the junction are eliminated by a tuning structure, and this is narrow-band. The 
solid lines in Figure g exhibit the frequency dependence of AG, a, e, and (p for the 25 MHz bands 
centered at four frequencies within the range commonly used with this feed. These particular 
data were derived from the source B0017-I-154; we obtained data for two additional sources, and 
the results agree well. The dashed lines are our adopted analytic expressions for the frequency 
dependence of the parameters. 

Near 1415 MHz the polarization is dual circular; the absence of a correlated cal means that 
the parameter ip has no meaning and, moreover, we cannot measure the position angle of linear 
polarization. The dependence of a on frequency is close to linear, which is what's expected for a 
turnstile junction. The variation of e is remarkably complicated, probably because of resonances 
in the tuning structure, and we do not have sufficient frequency coverage to characterize it. The 
scatter in 4> for the 1375 MHz spectrum simply reflects the uncertainty in the angle, which is 
large because e is small. Of course, AG simply reflects inaccurate relative cal values and not the 
properties of the turnstile itself. 

This work was supported in part by NSF grant 95-30590 to CH. 
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8. FIGURE CAPTIONS 

Fig. 1. — Polarizations of B1634+269 and B0017+154, derived from the calibration measurements 
for LBW. Position angles have not been rotated by Oastron- 

Fig. 2. — Mueher matrix parameters versus frequency for LBN, together with our adopted analytic 
approximations, ip is meaningless and not shown because this receiver does not have a correlated 
cal. 
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